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Various qualitative properties of solutions to the generalized Langevin 
equation (GLE) in a periodic or a confining potential are studied in this 
paper. We consider a class of quasi-Markovian GLEs, similar to the model 
that was introduced in ijj. Geometric ergodicity, a homogenization the- 
orem (invariance principle), short time asymptotics and the white noise 
limit are studied. Our proofs are based on a careful analysis of a hy- 
poelliptic operator which is the generator of an auxiliary Markov process. 
Systematic use of the recently developed theory of hypocoercivity [40] is 
made. 



1 Introduction 

In this paper we study various qualitative properties of solutions to the gener- 
alized Langevin equation (GLE) in M'* 



where V{q) is a smooth potential (confining or periodic), F(t) a mean zero 
stationary Gaussian process with autocorrelation function "f{t), in accordance 
to the fluctuation-dissipation theorem 



Here /3 stands for the inverse temperature and / for the identity matrix. The 
GLE equation ([!]), together with the fluctuation-dissipation theorem ^ ap- 
pears in various applications such as surface diffusion T and polymer dynam- 
ics [31] . It also serves as one of the standard models of nonequilibirum statistical 
mechanics, describing the dynamics of a "small" Hamiltonian system (the dis- 
tinguished particle) coupled to one or more heat baths which are modelled as 
linear wave equations with initial conditions which are distributed according to 
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{F{t)®F{s))^r^j{t-s)I. 
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appropriate Gibbs measures [SB]. In this class of models the coupling between 
the distinguished particle and the heat bath is taken to be linear and is governed 
by a coupling function p{x). The full Hamiltonian of the "particle+heat bath" 
model is 

H{q,p, (j), tt) = Hdp{p, q) + ^((^, vr) + j p{x)dq(j){x) dx (3) 

where H£)p{q,p) denotes the Hamiltonian of the distinguished particle whose 
position and momentum are denoted by q and p, respectively and 7^(0, tt) is the 
Hamiltonian density of the wave equation where^ and tt are the canonically 
conjugate field variables. The linear coupling in ([3]) is motivated by the dipole 
approximation from classical electrodynamics. By integrating out the heat bath 
variables and using our assumptions on the initial conditions we obtain ([!]), 
together with The memory kernel 7(i) in ([1} is given by the coupling 

function through the formula 

j{t) = Jimi^e'^'dk, (4) 

where p{k) denotes the Fourier transform of p{x) [20l |36] . 

The GLE has also attracted attention in recent years in the context of mode 
reduction and coarse-graininig for high dimensional dynamical systems |10j . 
One of the models that has been studied exte nsiv ely within the fr amework of 
mode elimination is the Kac-Zwanzig model [91 141] and its variants [26lfT5l[24l 
[3]. In this model, the heat bath is modelled as a finite dimensional system of 
N harmonic oscillators with random frequencies and random initial conditions 
distributed according to a Gibbs distribution at inverse temperature /3. The 
heat ba th c an be coupled either linearly or nonlinearly with the distinguished 
particle [5S] . Just as with model ^ , we can integrate out the heat bath variables 
explicitly. Passing then to the thermodynamic limit N — >■ +oo, we obtain 
the GLE The form of the memory kernel 'y{t) depends on the choice of 
the di stri bution of the spring constants of the harmonic oscillators in the heat 
bath jlOj . The Kac-Zwanzig model and its variants have proved to be very 
useful for testing various methodologies and techniques such as transition state 
theory ^il6j. 

The GLE ([T]) is a stochastic integrodifferential equation which is equiva- 
lent to the original infinite dimensional Hamiltonian system with random initial 
conditions. The infinite dimensionality of the original Hamiltonian dynamics 
with random initial conditions (or, equivalently, the non-Markovianity of the 
finite dimensional stochastic dynamics ([IJ) renders the analysis of this dynam- 
ical system very difficult. This problem was studied in detail by Jaksic and 
Pillet in a series of papers [221 EI] • In these works, existence and uniqueness 
of solutions as well as the ergodic properties of ([T]) were studied in detail. In 
particular, it was shown that the process {q,p = q} is mixing with respect to 
the measure 

Haidqdp) = ^g-^^^^^^'P) dqdp. 

To our knowledge, no information concerning the rate of convergence to equilib- 
rium for the non-Markovian dynamics ([l} is known for general classes of memory 
kernels. Ergodic theory for a quite general class of non-Markovian processes has 
been developed recently, see 12 and the references therein. 

A class of memory kernels for which more detailed information on the long 
time asymptotics of the GLE ([1]) can be obtained was considered by Eckmann, 
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Hairer, Pillct and Rey-Bellet in a series of papers [371 111 [3 [S]. Based on a gen- 
eralization of Doob's theorem on stationary, Markovian, Gaussian processes [1], 
it was observed in these works that when the memory kernel j{t) has a ratio- 
nal spectral density, then the GLE ([1]) is equivalent to a finite dimensional 
Markovian system. This system is obtained by adding a finite number of addi- 
tional degrees of freedom which account for the memory in the system. These 
auxiliary variables satisfy linear stochastic differential equations. As an example 
we mention the case where p{k) in ^ can be written as 



where p{k) = J2m=i Cm{—k)"^ is a polynomial with real coefficients and roots in 
the upper half plane. Then the Gaussian process with spectral density |io(fc)p 
is the solution of the SDE 

d , A dW 

where W(t) is a standard one dimensional Brownian motion. A related finite 
dimensional approximation of the infinite dimensional dynamics ([IJ has been 
introduced by Mori [30], see also [TT] and the references therein. Mori's tech- 
nique is based on a continued fraction expansion of the Laplace transform of 
the memory function j{t). 

Motivated by the above, in this paper we will consider finite dimensional 
approximations of the GLE. The general form of the Markovian approximation 
of (P) can be written as [24] [3 

Qm{t) = P,nit), Q,n(0)=g(0), (5a) 

Pm{t) = -VViQra{t))+X^z{t), P„(0)=p(0), (5b) 

z{t) = -P„,{t)X- Az{t) + CW, z(0) -AA(0,/), (5c) 

where z : R+ W" and A G R", A, C e M'"^". The fluctuation-dissipation 
theorem, which takes the form CC'^ — I3~^ {A + A'^), is assumed to be satisfied. 

In this paper we will consider ([5]) with A — (Ai, A2, . . . , A,„) and A diagonal 
with An = ai > 0. This amounts to approximating the memory kernel by a 
sum of exponentials, 

m 

^„(t) = ^A?e--l*l. (6) 

1=1 

It is expected that the results proved in this paper are also valid in the more 
general case ([S]). As remarked in [53], when A is invertible, the more standard 
Mori approximation |30| is equivalent to © after an appropriate orthogonal 
transformation. 

For this particular choice of A and A the SDEs ([S]) become (we drop the 



^To simplify the notation we consider {T} in one dimension. The results presented in 
this paper are valid in arbitrary finite dimensions. More details on the notation and on the 
multidimensional case can be found in Section |2. II and Remark 13.31 
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subscripts m for notational simplicity) 



Qit) = P{t), g(0) = <z(0), (7a) 

771 

-Vy(Q(i))+5]A,^(0, P(0)=p(0), (7b) 

-XjP{t) - ajz{t) + .y/2aj-^-Wj, Zj(0) - AA(0, /S"!) (7c) 
The process {Q(t), P{t)^ z{t)} is Markovian with generator 



i=i 

+ E (- W., - a,z, • V,^ + a,/3-i A,^ ) . (8) 

This is a degenerate second order elhptic differential operator of hypoelliptic 
type }19j . Convergence to equilibrium for models of the form ([7]) has been 
studied using functional analytic techniques [71 [5]. Similar results have also 
been proved using Markov chain techniques [281 137| . In this paper we present 
an alternative proof of exponentially fast convergence to equilibrium in relative 
entropy using the recently developed theory of hypocoercivity [40| . 

Several other results are also proved in this paper. We obtain sharp estimates 
on the derivatives of the Markov semigroup generated by —C. We prove a 
homogenization theorem (invariance principle) for ([7]) in a periodic potential 
and we obtain estimates on the diffusion coefficient. Finally, we study the white 
noise limit of the GLE ([T|), i.e. the limit as the noise F{t) in ([T]) converges to 
a white noise process. We show that in this limit the solution of (O converges 
strongly to the solution of the Langevin equation 

q^~V\q)--iq+^¥^W (9) 
and we obtain a formula for the friction coefficient 7 in terms of the coefficients 

The rest of the paper is organized as follows. In Section[5]we state our main 
results and we introduce the notation that we will be using. In Section|3]we prove 
exponentially fast convergence to equilibrium in relative entropy. In Section 2] 
we prove estimates on the derivatives of the Markov Semigroup generated by 
—C In Section [5] we prove the homogenization theorem. In Section [5] we study 
the white noise limit. For the reader's convenience, background material on 
the theory of hypocoercivity is summarized in Appendix |Xj Finally, the proof 
of geometric ergodicity of the process ([7]) using Markov chain techniques is 
presented in Appendix IbI 



P{t) = 

z,{t) = 

for j = l,...,m 
—C given by 

-C 



2 Statement of Main Results 

We will use the notation X := T"* x R'' x M''" and F M'' x M'' x R^™. 
We will also denote the process {q{t), p{t), z{t)} by x{t). When we study the 
dynamics Q in X the potential V{q) is periodic, whereas when x{t) E Y the 
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potential will be taken to be c onfining. The precise assumptions on the potential 
are given in Assumption 12.11 below . 

Our fist result concerns the ergodicity of the SDE ([7]) in X or in y. To prove 
the ergodicity of the SDE in Y we need to make the following assumptions on 
the potential. 

Assumption 2.1. 

(i) V{q) G C^(R'') is a confining potential. 

(ii) {WqV,q) > (yV{q) + /3||<z|P for some suitable /3,a > where (•,•) and || • || 
denote the Euclidean inner product and norm, respectively. 

(Hi) There exists a constant c such that ||V^y|| < c, where \\ ■ \\ denotes the 
Frohenious- Perron matrix norm. 

The invariant measure fip{dqdpdz) — p(p, q,z) dpdqdz of the process ([7]), 
whose density satisfies the stationary Fokker-Planck equation C* p = 0, is known: 

p(g,p,z) = ie-'^(^l^'l'+^(«)+^ll^ll'), (10) 

where 2, is the normalization constant. This invariant measure is unique and 
the law of the process ([7]) converges exponentially fast to pfi. 

Theorem 2.1 (Ergodicity). The process with x{t) e X and V{q) € C^iT^) 
is geometrically ergodic. The same holds true when x(i) S Y , provided that the 
potential V{q) G C'^(M'^) satisfies Assumption \2. 1\ 

The proof of this theorem, which is based o n M arkov cha in-t ype arguments 
and which is similar to the proof presented in [37], see also [2S], can be found 
in Appendix IB] 

We can prove exponentially fast convergence to equilibrium using tools from 
the theory of hypocoercivity |3D]. We will use the notation K. := Ker(£) and 
77^ for the weighted Sobolev space with respect to pp on either X or Y . 

Theorem 2.2. Let —C he the generator of the process x(<) G X , the solution 
of (I7|). Then there exist constants C, A > such that 

The sam e ho lds true when x(t) G Y, provided that the potential satisfies As- 
sumvtion \2. lY i) and (Hi). 

Using the tools from ['TO] we can prove exponentially fast convergence to 
equilibrium in relative entropy. The relative entropy (or KuUback Information) 
between two probability measures /i and v with smooth densities / and p, re- 
spectively, is defined as 

Hp{f)= y/iog(^Q dx. 

We will measure the distance between the law of the process x(f) at time t and 
the equilibrium distribution in relative entropy. 
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Theorem 2.3 (Convergence to Equilibrium). Let ft be the law of the process 
x(t) at time t and assume that Hp(fo) < +00 and V{q) £ C^(T''). Then there 
exist constants C, a > such that 

Hp{ft) < Ce-"*HM. 

The same holds true when x(t) € Y , assuming that Hp{fo) < +00 and provided 
that the potential V{q) satisfies Assumvtion \2. lY i) and (Hi). 

Remark 2.1. In view of the Kullback inequality 

\\\ft^p\\L-<Hp{ft), (11) 

Theorem implies that, for initial data with finite relative entropy, we have 
exponentially fast convergence to equilibrium in . 

The proofs of Theorems 12.21 and 12.31 are presented in Section [3l 
Estimates on the Markov semigroup and its derivatives associated to the 
Lan gevin equation can be proved using an appropriate Lyapunov function jl8l 
113) . In this paper we use similar techniques to obtain estimates on the Markov 
semigroup and its derivatives for the generalized Langevin equation, equa- 
tion 0. We introduce Cfc, fc = 0,1,2 with Co = A, Ci = [A,B] and C2 = 
[Ci, B]. We will use the notation := i^( • ; iip{dx.)) where • is either X or Y. 

Theorem 2.4 (Estimates on Derivatives of the Markov Semigroup). Let —£ be 

the generator of the process x(t) g X , the solution of ^ with V{q) G C^(T'^). 
Then the Markov semigroup e^*^ satisfies the bounds 

\\Cke-'''\\Li^Ll<C^, k ^0,1,2 and te [0,1]. (12) 
t 2 

The same holds true when x(t) G Y , provided that the potential V{q) satisfies 
Assumvtion \2. lY i) and (Hi). 

Remark 2.2. This result can also be obtained by applying Theorem lA.Sl Malli- 
avin calculus-type arguments show that estimate ilS^) is sharp. 

When the potential V{q) is periodic, the particle position, appropriately 
rescaled, converges weakly to a Brownian motion with a diffusion coefficient 
which can be calculated in terms of the solution of an appropriate Poisson 
equation. Results of this form h ave been known for a long time for the Smolu- 
chowski (overdamped) equation [SH Ch. 13] as well as for the Langevin dynam- 
ics |33[ 114] . In this paper we prove a similar result for the generalized Langevin 
equation. We will use the notation (ff- := (p ■ e, p'^ := p ■ e, where e denotes an 
arbitrary unit vector in M.'^. 

Theorem 2.5 (Homogenization). Let x(t) G X be the solution of ^ with 
V{q) G C°°{T'^) with stationary initial conditions. Then the rescaled process 
ql{t) := e- eq{t/e^) converges weakly on C([0,T],R) to a Brownian motion with 
diffusion coefficient D with 

m 
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where (jf £ L'j^ is the unique, smooth, mean zero, periodic in q solution of the 
Poisson equation 

= (14) 

on X . Furthermore, the following estimates hold 

i—l 

Let q{t) be the solution of the Langevin equation © and let q'^{t) := q{'yt). 
It is well known that this rescaled process converges in the overdamped limit 
7 — > +00 to the solution of the Smoluchowski equation [3TJ Ch. 10] 



q = -VV{q) + V2/FiPF. (16) 

Similar results have also been proved in infinite dimensions |38| . In this paper 
we prove a similar result for the convergence of solutions to the GLE to the 
Langevin equation in the strong topology and obtain a formula for the friction 
coefficient that appears in the Langevin equation. 
Consider ([T]) with the rescaled noise process 

F^it) := -^/(t/e), (17) 

which is a mean zero stationary Gaussian process with autocorrelation function 

Y{t) = -^^(t/e). (18) 

For the memory kernel (H]), 7'^(i) becomes 

rn A, 2 

7^(0=^X^^1*1. (19) 

Consequently, the rescaled noise process (|T7| is obtained by rescaling the coef- 
ficients in ([7]) according to \j ^ aj ^ Under this rescaling the SDEs 
become 

r m = p{t), 

Pit) = -d,V{q)+j:t,^^z,it), (20) 



[m = -^Pt - f ^.(0 + V^^^- z = l,...,m. 

Theorem 2.6 (The White Noise Limit). Let {q{t),p{t),z{t)} e X be the solu- 
tion of 1120]) with V{q) G C^iT^) and initial conditions having finite moments of 
any order. Then {q{t) , p{t)} converge strongly to the solution of the Langevin 
equation 

Q{t) = P(t), 

Pit) = -^,v{Q,)~j:ti(im~\f^w^ (2^) 
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Consequently, the process {q,p} converges weakly to the solution of the Langevin 
equation 

-Q{t) = Pit), (22) 



P{t) = -VgV{Qt)-jP+^/2^W, 
where the friction coefficient 7 is given by the formula 

a,: 



7 = E' 



2.1 Notation 



For x{t) = (g,p,z) £ Y := R'^ x M.'^ x M*" or x(t) e X := T'- x 
consider the operator C: 



m 

with kernel /C := KerC. The density of the invariant measure of the process 
x(i) is 

(24) 

where | • | denotes either the Euchdean or the matrix norm. In (P^ . V is the 
gradient (or the derivative when d = 1) and A the Laplacian. denotes the 
Hessian and if O is an operator then O* is its adjoint in := ■ ; /i^((ix)). 
Furthermore, we wih use the notation L = —C so that L is the actual generator 
of the process. 
Define 

m 

B = -pVq + \7qV\7p - J2 ^3 (^J^P " P^-, ) • (25) 

We easily check that B* = —B. If m = 1 then Ai = —d^i (derivative with 
respect to the i — th component of z) so that A* — —Zi + dz^ and we can write 

d 

/: = B + ^A*Ai =B + A*v4, (26) 
1=1 



where A is intended to be the row vector of operators {Ai, ...,Ad) ( the same 



for A*). More precisely, if m = 1 then: A : — > <^ R'^ , B : 



Ll, [A*, A] : Ll — ^ L% being [A*, A] := J2j=M*^Aj]; on the other hand 
[A, A*] : L^p — > L^p^R!^ ® K.'' is a matrix whose «j-th component is given by 
[A, := [Ai, A*]; in an analogous way [A, A] : — > (g) M'' (g) E** is a 
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matrix with [A,A],j [A,, A^]; finahy C := [A,B], C : Lf, — ^ ^ is a 
vector of operators, Ci = [Ai, B], i = l...d, and the same holds for C2 := [C, B], 

If m > 1 then (PBl) becomes 

m d 

£ = 5 + 5]^ a:,. A,, (27) 

with = — 9z;. i.e. derivative with respect to the j-th component of Zi. We 
will use the notation 

C = B + A*A, (28) 

meaning either or (|27p . 

As for the norms, unless otherwise specified, || • || indicates the norm of L^, 
||-||i = ||A-|p + ||C-|p + ||C2-|pisa sort of homogeneous H^{Y\ ^/3(dx)) norm 
and II • 11^1 = II • 11^ + \\A ■ Ip + ||C • |p + IIC2 • |p is the usual inhomogeneous 
one. The inner products in these Hilbert spaces are denoted by (•,•), (•,•)! ^'iid 
(•, •)//!, respectively. 

Remark 2.3. It is a general result that if C is accretive and C is maximally 
accretiv e th en C is the generator of a contraction semigroup, so e~^^ is well 
defined {I^ Ch. 5]. 

L is accretive since, as already noticed, (Cf, f) = ||A/|p > V/ in the domain 
of C To prove that C is maximally accretive we show that the operator C* + vl 
is infective for all v > Q. Since the operator is linear, we just need to show that 
Ker{C* + vl) = {0} for all v > Q. So suppose f £ Ker(^C* + vI), then 

C*f + iyf = -Bf + A*Af + i.f = 0. 

Consequently 

i-BfJ) + \\Afr + u\\fr = 0. 

Therefore 

0<\\Aff = -i^\\ff^f = 0. 

Notice that when v = 1, from the above we conclude that A* A is self-adjoint and 
not only symmetric \S5[ Thm VIII. 3]. From now on by C we mean its closure. 

3 Convergence to Equilibrium 

3.1 Hypocoercivity 

Background material on hypocoercivity is presented in Appendix [Xj 

Definition 3.1 (Hypocoercivity). With the same notation as in Definition \A.l\ 

and assuming that the operator T generates a continuous semigroup, such an 
operator is said to be X-hypocoercive on % if there exists a constant k > such 
that 

II ||^< Ke"-^* II h 11^ V/i e ?i andt>0. 
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We say that an unbounded hnear operator S' on is relatively bounded with 
respect to the (hnear unbounded) operators Ti,...,Tn if 'DiS) C (n2?(T,)) and 
3 a constant a > s.t. 

yh e V{S), \\Sh\\ < a{\\Tih\\ + ... + \\T„h\\). 

The basic idea employed in the the theorems that we are going to use is to 
appropriately construct a scalar product on by adding lower order terms 
and then use the fact that hypocoercivity is invariant with respect to a change 
of equivalent norms, whereas coercivity does not enjoy such invariance. Finally, 
notice that S, the class of Schwartz functions is dense in L^, hence it is dense 

in D{A) n D{B). This guarantees that all the operations performed on these 
(unbounded) operators are well defined. 

Set m = l = d, a = X = P = l. The first two commutators are 

Ci = C = [A, B] = dp and C2 = [C, B] = - dg. (29) 

Hence the operator is hypoelliptic [19) . Furthermore, 

[A, A]^0 [A, C] = [A, C2] = 0, 
[A,A*]^Id [C,A*] = [C2,A*] = -Id, 

[C, C*]^Id [C; ,C2]^-Id~d^ V, (30) 
where Id is the identity operator. 

3 . 2 Proof of Theorem 

Proof. We will use Theorem IA.2I . To this end, set 

P ^ A*A + C*C + C2C2 

and notice that Ker{P) = K, =: KerC contains only constants; in fact 

Ker(P) = Ker{A* A)r\Ker{C* C)<r\Ker{C*2C2) = Ker{A)nKer{C)nKer{C2). 

To show that K. = Ker{A*A) n Ker{C*C) n Ker{C^C2): the inclusion D is 
obvious. For the other inclusion: if /i S /C then || + ||C/ip + \\C2h\\^ = ^ 
Ah = Ch = C2/1 = 0. 

The above mentioned theorem re quire s two sets of hypotheses to be fulfilled. 
Hypothesis 1,2 and 3. in Theorem I A . 2 1 are quantitative assumptions, which are 
satisfied in our case with N = 2, Co — A, Ci — C, Ri — R 2 ^ 0, R3 — [C2, B] 
(this is to have C3 = 0) and thanks to Assumption I2.1f iii'). Hypothesis 4. 
requires, in our case, for the operator P to be K-coercive on /C^ = L^p/JC. The 
coercivity of P is equivalent to 

\\Ahf + \\Chf + \\C2hf > K\\hf, 

that is, more explicitly, 

WV.hf + \\V.phf + II (V, - V,) hf > K\\hf. 
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Using the fact that |la - 6|p > we have 

II V,/if + \\Wph\\^ + II (V, - Wg) hf>^ {W^zhf + + WWghf) 

so we just need 

||V,/i||2 + ||Vp/i||2 + ||V,/i||2>K||/if 
to hold true. Since /i^ is a product measure, we only need to verify that 



holds true for some constant /i, where the notation (h) := J he ^^i^'^i has been 

used. It is a standard result that if V{q) e C^{R'^) is such that e^^^'^/Z is a 
probability density and 

^^-Ay(,)in-+oc (31) 



then e~^^^V Z s atisfies a Poincare inequality (see, e.g., [40l Thm. A.l] ). From 
Assumption I2.ir iii). Condition ([3T|) is satisfied. We can conclude that there 
exist a scalar product ((•,•)) inducing a norm equivalent to the inhomogeneous 
norm of and a constant A > such that C is coercive in this norm: 

Vh e L'l/IC, {{h, Ch)) > X{{h, h)). 

This implies that £ is hypocoercive in this norm, hence it is hypocoercive on 
L'^/JC endowed with the || • ||/fi norm: 

\\e-'''-'ho\\m <Ce-^'\\ho\\m. (32) 

□ 

Remark 3.1. The orthogonal space to JC is the same with respect to both the 
and the norms; moreover, since P is coercive, these two norms are 

equivalent. 

Remark 3.2. Theorem \A.3\ in Appendix A allows us to state a similar result 
when the initial datum is in L^. In fact 

||e"*^/i|ki < ^||/i||, te{0,l]. (33) 

t 2 

So, putting together i32\) and i33\} we get, for < tg < t, tQ < I: 
||e-*^/io|ki =||e-(*-*»)^e-*«^/.o||Hi - \\e-^'~'o)c^^j^^ 

< ce-^(*-*«)||/jiJ|^i < ce-^(*-*°)||e-*°^/io||i/i 

g-A(t-io) 

< C 5 ll^oll- 

'-a 

Remark 3.3. The proof is identical when m,d> 1. In this case we can think 
of A as a matrix of operators, see \21^ . 
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3.3 Proof of Theorem [273] 

Proof. Let ft denote the law of the process x(i). We set ft = pht- Then ht 
satisfies the equation 

dtht = Bht - A*Aht . (34) 
We apply Theorem lA. 41 to the operator F = —B + A* A with 

A = -d,, Ci = -dp, C2 = Z2 - Id, i?2 - -5,. 

Furthermore Asuumption l2.1l fi) and (iii) together with the HoUey-Strook per- 
turbation Lemma imply that Z~'^e~^^'^'' satisfies a Logarithmic Sobolev Inequal- 
ity (LSI). 

Hypotheses 1,2 and 4 are automatically satisfied. We put C2 = dq and we 
added the remainder R2 in order to fulfill hypothesis 4. Hypothesis 3. is satisfied 
on account of Assumption 12. II (iii). Now consider the relative entropy Hp{f), 

Hp{f) = J flog (^^^ dqdpdr ^ Jhloghdp f^ph (35) 

and the Fisher information Ip{f) 

IpU)= J f \^\og{h)\^ dqdpdr = J h\V\ogh\^dp f^ph. (36) 

Then if the initial datum has finite relative entropy, we obtain that 

Hpift) = 0(e-*") (37) 

for some a > and for t > 0. If the initial datum has also finite Fisher 
information then 

lp{ft)^0{e-"^), (38) 

as well. □ 

Remark 3.4. We remark that S38\} . together with the LSI, implies (5' 



Rem ark 3.5. In viev of the LSI, it is interesting to notice that, applying The- 
orem \A.5l we get the following bounds 



c 



h{t)\Cklogh{t)\'dp< J hologhodp, (39) 

for k — 0,1,2 and c an explicitly computable positive constant. 

4 Bounds on the derivatives of the Markov semi- 
group 

Throughout this section we well be using the notation u = c^^'-'Uq. We introduce 
the Lyapunov function 

F{t) = a^tWAuf + ait'^WCuf + a2t^\\C2u\\'' 

+bat^{Au, Cu) + t%i{Cu, C2U) + b2\\uf, t € (0, 1], (40) 
where aj, bj, j = 0,1,2 are positive constants to be chosen. 
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Lemma 4.1. There exist constants aj,bj,j = 0, 1,2 such that the time deriva- 
tive dtF of the Lyapunov function along the semigroup is negative. 

Proof. We will calculate the time derivative of each term in (|40)) separately and 
using the explicit relations ([30]) : 

dt\\uf = -2{Cu,u) = -2\\Au\\'^, 

dt{Au,Au) = -2{Cu,Au)~2\\A*Auf ^~2{Cu,Au)-2\\Au\\^ ^2\\A^u\\'^, 

dt{Cu,Cu) = -2\\ACu\\^ -2{C2U,Cu), 

dt{C2u,C2u) = {{2 + dlV)C2u,Cu)-2\\AC2uf + 2{Au,C2u), 

dtiAu,Cu) = -2{A'^u,ACu)-{Au,Cu)~\\Cu\\^ ~{Au,C2u), 

dt{Cu,C2u) = -\\C2u\\^-2{ACu,AC2u) + 2\\Cuf + {Cu,Au). 

Putting everything together we obtain 

dtF{t) = -2aot\\A'^u\\^ ^2ait^\\ACu\\'^ -2a2t'^\\AC2u\\'^ 
-2bot'^{A^u,ACu) - 2bit'^{ACu,AC2u) 
+{-2aot + ao- 2b2)\\Auf + (Soii^ + 2bit^ - bot'^)\\Cuf 
+ {5a2t^ - bit%C2uf + (26ot - 2aat - bot^ + bit'^){Au, Cu) 
+ (461^3 - 2ait^ + 2a2t^){Cu, Csu) + (2a2i^ - bot^){Au, Caw). 

Now we estimate the sum of the first and of the second line (i.e. the sum of all 
the terms where A^, AC and AC2 appear). For t E (0, 1] we have 

(mil + (1311)2 < ~2aot\\A^uf + 2bot^\\A'^u\\\\ACu\\ 

+2bit^\\ACu\\\\AC2u\\ - 2ait^\\ACuf - 2a2t^\\AC2uf 
< ^2aot\\A^uf + blt\\A^uf + t^\\ACuf - 2ait^\\ACuf 
+blt^\\ACuf + t^ll - 2a2t'^\\AC2uf. 
Similarly for the sum of the remaining terms (those with A, C and C2)we have 

(SIDs + (I3ID4 + (HIDs < {-2aat + ao - 2b2)\\Au\\^ 

+{2bot + 2aat + bat'^ +bit'^)\\Au\\\\Cu\\ + {2a2t^ + bot^)\\Au\\\\C2u\\ 
+{3ait^ + 2bit^ - bat^)\\Cuf + {5a2t^ - bit^)\\C2uf 
+ {Abit^ + 2ait^ + 2a2t^)\\Cu\\\\C2u\\ 
< {-2aat + aa- 2b2)\\Auf + al\\Auf + \\Cuf 

+alt'^\\Au\\^ + t^\\C2u\\^ + jbl\\Au\\^ + ^^\\C2uf 
+2blt^\\Cu\\'^ +t^\\C2u\\^ + aft^WCuW^ +t^\\C2u\\^ 
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alt''\\Cuf + t'>\\C2uf. 



Choosing the constants m such a way that 62 ^ oo ^ ^ ^2 > 1/c 



where c is a constant depending on the bound on the second derivative of the 
potential, we obtain that dtF < £ (0, 1]. □ 



Proof of Theorem \2.4\ We use the previous Lemma to deduce 

aot\\Au\\^ + ait^WCuf + asi^UCsuH^ 
+bot'{Au,Cu)+t%i{Cu,C2u) + b2\\uf < 62IIU0IP. 

This, in turn, imphes that 

\\VMf = \\Auf<^\\uo\\^ 

\\v,ur = \\cur<^\\uor, 



2 



K 



=^ ||V,U||2 < Jlluoll', 

where k is an exphcitly computable positive constant. The previous inequahties 
are justified by the fact that 

aot\\Auf + ait^WCuf + a2t^\\C2uf + bat^{Au, Cu) + t%i{Cu, C2U) 

and the second line is positive thanks to the choice of the constants we made. □ 
Remark [33] holds also in this case. 



Remark 4.1. From the estimates similar estimates onA*e~*^ , e~^^ A*, 
C*e~*^ , e~*^ C , 6*2 e^*^ and e~*^ C* follow, where* and* stand for either 
the L^-adjoint or nothing. In fact: 

(z) {Ae-'^f,g) = {f,e-'^'A*g) < || Ae-*^/|l II5II < ^ll/llbH 

(/, e^*^'v4*,g) < -^11/11 ||g||, choose f — e^*^'A*g and the result on 
e~*^* A* follows. 

(ii)Vsing [A,A*]^Id we have P*e-*^Mof = \\A*uf = \\Au\\^+\\u\\^, hence 
the estimate for A*e~*^. Taking the adjoint as in (i) we get the result for 

(Hi) For Ae"^^ we can just repeat the proof we wrote for Ae^*"^ , since the only 
thing that changes when considering C* is the sign of B, which doesn't 
play any role in the proof. Now, by acting as in (i) and (ii), we obtain the 
results for e^*'~ A* ,A* e^*'-' and e~*^A. 
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5 The Homogenization Theorem 



In this section we prove Theorem 12.51 The proof of this theorem is based on 
standard techniques, namely the central hmit theorem for additive functionals 
of Markov processes [531I171E3], which in turn is based on the martingale central 
hmit theorem [HI Thm. 7.1.4]. In order to apply these techniques we need to 
study the Poisson equation 

l^u = f. (41) 
The boundary conditions for (j41l) are that u ^ and it is periodic in q. 

Proposition 5.1. Let f £ L^nC°°(X) with j-^ f iij3{(hC). Then the Poisson 
equation (|4ip has a unique smooth mean zero solution u £ L'^f] C°°{X). 

The proof of Theorem 12.51 follows now from the above proposition. 

Proof of Theorem \2.5\ To simplify the notation we present the proof for d = 1. 
When d> 1 the same proof applies to the one-dimensional projections q'^ — q-e. 
In this case the diffusion coefficient D is replaced by the projections of the dif- 
fusion tensor D'^ = De ■ e. 

We consider the process x(i) on X with stationary initial conditions. For 
non-stationary initial conditions we need to combine the analysi s pr esented be- 
low with the exponential convergence to equilibrium, Theorem 12.21 Since p € 
L^nC°°(X) and centered with respect to the invariant measure ij,p{dx.), Propo- 
sition [O] applies and there exists a unique mean zero solution (j> £ L'j^D C°°{X) 
to the problem 

£0=p. (42) 

We use Ito's formula to obtain 



d(j) = C(f)dt + Y, \l'iajl3-^d,^(j3dWj. 

We combine this, together with and the equations of motion to deduce 
q,{t) eqit/e") 

= eq{0) + e / p{s) ds 
Jo 

= eq{0) ~ e{^{q{t/e^),p{t/e^),z{t/e^) - 0(g(O),p(O), z(O)) 

+eE / j2a,p-^d.^4>dW,{s) 
=: eW + M\ 

Our stationarity assumption, together with the fact that (p £ L^, imply that 

\\R'\\ < C. 

To study the martingale term NP we use the martingale central limit theorem [51 
Thm. 7.1.4] or ;34, Thm 3.33]. We have that Af^(O) = 0, that M\t) has 
continuous sample paths and, by stationarity, that it has stationary increments. 
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Furthermore, by stationarity and the fact that the Brownian motions Wi{t), i 
1, . . . m are independent, we deduce that 



lim(M,^) =2^a,/3-'l|5..0|pi in 
1=1 



The above calculations imply that the rescaled one dimensional projection qdt) :— 
eq{t/t^) converges weakly in C([0, t]; M) to a Brownian motion ^2DW{t) where 

m 

D = 2^a,/3-i||a,,0f . (43) 
1=1 

Remark 5.1. Notice that when d > 1 the convergence of the one dimensional 
projections does not imply the convergence of the process qe{t) = eq{t/e^). The 
proof of this result, which is also based on the analysis of the Poisson equation, 
is very similar and it is omitted. 



To prove estimate (|T5)) , we first show the upper bound and than the fact that 
the diffusion coefficient is bounded away from zero (when we consider periodic 
solutions only). We set (p — gi + -^Zi and use the Poisson equation (|^ to 
obtain 

from which we obtain the estimate 

m 

a.p-'\\dzM^ < J2a,|3-'\\^.,g^\\'^{£g^,g^) 



Consequently < j-. From this we obtain the following estimate on the 

diffusion coefficient D 

m m 

1=1 ' i=l * 

2 ™ / 1 

A ^ 

< !^ 

1=1 ' 

The fact that _D > is easily seen by contradiction. If = then, by (|33]), 
ll^^^f/ilP = Vi. Hence — 4>{q,p) and 

m 

C(/) = -pdqcj) + dqVdpcj) + ^ XiZidp4> = p- 

i=l 
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Multiplying both sides by 1"^ and then integrating with respect to Zi we get 

+ y Xizfe^'^/'^dz^ + ^ y XtZiZjdp(j)e''^/'^dzi 
= y pz e^'^ ^'^dzi, 

from which we conclude that Xidpcj) = for all i. Hence (j) = (t>{q)- By the same 
reasoning we get that —pdqcf) = p, which does not have a periodic solution. □ 

Proof of Proposition \5.1\ 

The scheme of the proof, which is similar to the proof of [531 Lemma 2.1], is 
as follows. We consider the Poisson equation £0 — f where f € LpH C°°(X) 
and centered with respect to the invariant measure fip{dx). 

1. We consider the modified problem 

Xu^ + Luy + vPu^ = / A, > (44) 

where P is a regularizing operator (i.e. the operator £ + P is uniformly 
elliptic) so that the Lax-Milgram theorem applies; hence the weak solution 
to (|44l) is unique Vi^ > and letting — 0, by the uniqueness of the weak 
limit, we get existence and uniqueness of the solution to 

Aw + = / A > (45) 

(notice that — £ is the generator of a Markov semigroup so by Hille-Yosida 
theorem the set {A £ P : A > 0} is contained in the resolvent of the oper- 
ator — £) 

2. Set ll,\u = Am + Lu ~ Xu + f (the last equality has to hold in distribution 
if we want u to satisfy (1411) ) 

^u^ Cl\Xu + /) - XC^^ + Cl^f 
and defining C^^ f = h we get (jld — C^^)u — h, h — h/X. 

3. Once proven that is compact we are done; in fact in this case the 
Fredholm Theorem applies so either the solution to {■jid — C~^^)u = h 
exists and is unique (and hence, by construction the solution to (j4ip is 
unique) or (^/d— £^^)u — admits a nonzero solution. We can rule out 
the latter option because {\ld — jC-'^^)u = 4^ Lu = 0; since we know 
that KerC contains only constants and we are asking for u to have mean 
zero we can conclude that Cu = Q <^ u = Q and we are done. 

Now the details. Choose P — C* C + G* G where G = Vg and notice that now 
B = G*G-C*G+A*G-G*A. Let us check that the hypothesis of Lax-Milgram 
Theorem hold: 

A(u„ u,) + \\Au4^ + iy\\Cu,\\^ + iy\\Gu4^ 
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and 

\{uu,v) + {Buu, v) + (j4ii,y, Av) + v{Cuu, Cv) + v{Guu, Gv) 

< AK,«) + ||Cu,||||Gt;|| + ||GM,||||Ct;|| + 
+i^\\Cu,\\\\Cv\\ + i^\\Gu,\\\\Gv\\ + \\Gu,\\\\Av\\ + \\Au,\\\\Gv\\ 
< a\\u,y\\m\\v\\Hi, 

with a a positive constant depending on ly and A. Hence the weak solution to 
(|44)) exists and is unique. Moreover, from (|44|) we also have that 



A||u,||2 + pu,||2 + j.(||C«,||2 + ||G^.,||2) = < ||/||||u,||, 

hence A||ui/|| < c, < c and i/||Cui/||^, j^||Gui/|p < c. Now taking as 

test function we get that < c||/|| and we can let — > obtaining that 

Ui, converges weakly to the solution of ([^5)1 . Recalling that u — C^^f and that 
u the weak limit of {ui/}, we can write 

WJ^x'fWm = \\u\\m < liniinf \\u,\\m < c\\f\\ 

and using the fact that Hj^ is compactly embedded in we can conclude that 
the resolvent is compact. This completes the proof. □ 



6 The White Noise Limit 

Throughout this section C denotes a generic constant, which is independent of 
e. To simplify the notation we present the proof in one dimension. The proof 
is exactly the same in arbitrary dimensions. Let {Q{t), P{t)) g T x M be the 
solution to the system (|2T|) . Then 



I q{t) ~ Q{t) \<\ q{0) - Q(0) I + /* I Pis) ~ P{s) I ds. 

Jo 

Setting 6i = A| /ai we get that 

m m , 

p{t) - P{t) = -d,v{q) + d,v{Q) + J2 UP{t) - Pit)) - V^E ^^»(*)- 

4=1 1=1 * 

Hence, for any r > 2, 

77(T) := E sup {\q{t)-Q{t)r + \p{t)-P{t)n 
te[o,T] 

< Cr-^ [ E sup \ qis) - Q{s) r dt 
Jo se[o,t] 

+G If^eA T^-' f E sup I p{s) - Pis) r dt 
Vj=i / "'0 selo,t] 

™ /A A"" 

+Ge5 V - E snp \ z^it) - z,(0) T \ p(0) - P(0) T +E \ g(0) - Q(0) T 

^ t6[0,T] 



i=l 
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From Gronwall's Lemma we get 



viT) < C 



V p E snp |z,(t)-z,(0) 
^ ^ ^ I te[Q,T] 



£26 



CT 



and the result now follows from Proposition 16. II 

Proposition 6.1. With the same notation and assumptions of Theorem \2.6[ the 
following estimate holds true: 

^(T) sup |z,(Or+-Bsup \ p{t) \-< C ( 1 + T + -T + J -T] 

.^-^ te[o,T] te[o,T] \ e V e / 



+C 1 



A a + y/a 



/ mdt+E\p{Q)y +y^E\zm 

■^0 ^=l 



In particular 



y^E sup \z^{t)\'-<Ce-\T + VT)e 



£; sup I - z,(0) r — ^0 Vr>2. 

t6[0,T] 



Proof. We use Ito's formula to deduce that 



1 ft ^ rn „t 

-pW- {d,V{q{s))p{sY-')ds+ — Y, Kp{sr^'z,{s)ds 

Jo V £ j^;^ -'0 

-z,{OY~- f z,{srds-^ f\p{s)z,{sY-^)ds 



+ - f z,{sY-^ds+\l^ f z,[sY-^ dW,[s). 

Setting A = maxjA^, i = l...m}, using the boundedness of dqV on the torus and 
Holder's inequality, we have 



E sup \p{t)\''<CT + C (l + [ E sup 



+ ^mA /" V £; sup I z,{s) r dt + £: I p(0) 



19 



Similarly, 

V-B sup I r < C'^A / sup \p{s)\'^dt 
~{ te[o,T] ye Jo se[o,t] 

X /"9 /"^ ™ 

+(7(4 +2- + J—) y sup I zi{s) r 



i— 1 



In the above we used the following estimate, which is a consequence of the 
Burkholder-Davies-Gundy inequality 



E sup 

te[o,T] 



Zi{sy-UWi{s) 



< E 



r I zi{s) i^^'-i) ds 

Jo 

zi{s) 1^'- +i)ds 

Jq 



< CE 



< CE sup I Zi{t) r + Cri 
te[o,T] 



< CT^ / E sup I Zi{s) r dt + CT^. 
Jo se[o,t] 

Putting everything together we obtain 

m 

y^E sup \zi{t) r+£; sup \p{t) r 

~t te[o,T] te[o,T] 



< C <T + m-T + m^ 



+ C{1 + ^) [ E sup \p{s) \^ dt 

Jo se[o,t] 

+ C( A + ^+ J^) rf2E sup I Zi{s) r dt 
ye e V e Jo ^ s6[o,t] 



j=l «6[0,t] 

m 

+ i;|p(o) r+^£;h,(o) r 



i=l 



so that 



Jo 

Applying Gronwall's Lemma gives then 

e^^iT) < C{T + yT)ei-ie^^^'"' 
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and, since £2 J^^i ^^^Pt£[Q t] I ^ii^) T— we get the result. 

□ 

Acknowledgments. The authors are grateful to Dr Karel Pravda-Starov for 
many useful comments and for an extremely careful reading of an earlier version 
of the paper. 



A Hypocoercivity 

In this appendix we recall some of the main results from the theory of hypocoer- 
civity, as presented in 140 . Throughout this Appendix we will use the notation 
introduced in Section l2.ll 

Definition A.l (Coercivity). Let T be an unbounded operator on a Hilbert 
space % with kernel /C. Let % be another Hilbert space continuously and densely 
embedded in JC^ . The operator T is said to be \-coercive on T-L if 

{Th,h)^>X\\h\\^^ yhelC^nD{T). 

The following Proposition gives an equivalent Definition of coercivity. 

Proposition A.l. With the same notation as in Definition \A.l\ T is X-coercive 
on T-L iff 

II er'^^h \\^< e"^* II h 11^ V/i e H andt> 0. 

Theorem A. 2. Let L be an operator of the form C — A*A + B, with B* = —B, 
K — KerC and assume there exists N such that 

[C,-i,B]=C,+R, l<]<N+\, Co^A,Cn+i^O. (46) 

Consider the following assumptions: for k = Q, N + \ 

1. [A, Ck] is relatively bounded with respect to {Cj}o<j<fe '^'^'^ i^i^}o<j<k-i 

2. [C'k,A*] is relatively bounded with respect I and {C'j}o<j<fc (here I indi- 
cates the identity operator on L'j^) 

3. Rk is relatively bounded with respect to {Cj}o<j<fe-i '^^^ {^J^}o<j<fc-i 
4- '^jLo K-coercive for some k > 

If Assumptions 1 — 3 are satisfied then there exists a scalar product ((•, •)) on 
defining a norm equivalent to the usual norm and such that 

N 

WheH^/K., {{h,£h))> Kj2\\(^jhf (47) 

3=0 

for some constant K > 0. Furthermore, if Assumption 4 is satisfied, then there 
exists a constant A > such that 

Vh e H^/JC, {{h, Ch)) > X{{h, h)). 
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In particular, C is hypocoercive in H^/JC, i.e. 

\\e-*^\\H^/K^H^/K<Ce-^' 

for some C, A > 0. 

Theorem A. 3. With the same notation as in Theorem \A.^ if Assumptions 
1-3 are satisfied then 

||Cfee-*^/i|l < Vfc = 0,...,iV 

for all functions h £ L^. 

Theorem A.4. Let V{x) e C^iR"^) such that n{dx) = e^^'-'^Ux is a probability 
measure on R'^ and assume that C generates a semigroup on a suitable space of 
positive functions. Let {Aj}l<^j<^]^J and B be first order differential operators 
with smooth coefficients, with B = —B*. Assume there exists N £ N such that 

[Cj^i,B]=Cj+Rj l<j<N + l, Co = A,Cn+i^O. 

If, for < k < N + I the following assumptions are fulfilled 

1. [A^Ck] is pointwise bounded with respect to A. 

2. [Ck,A*] is pointwise bounded with respect to I and {Cj}o<j<fc' 

3. Rk is pointwise bounded with respect to {Cj}^^j^i^_^. 

4. [A,Ck]* is pointwise bounded relatively to I and A. 

5. there is a positive constant A > such that X^fc ^fc^fc — pointwise on 

(I is the identity matrix on W^). 

6. The probability measure fi satisfies a logaritmic Sobolev inequality. 

Then the Kullback information and the Fisher information decay exponentially 
fast to zero. 

Theorem A. 5. With the same notation as in Theorem \A.4\ let V{x) S C^(]R'^) 
be such that fi{dx) = e~^^'^^dx is a probability measure on R" and assume that C 
generates a semigr oup o n a suitable space of positive functions. If Assumptions 
1 — 4 0/ Theorem \A.4\ are fulfilled, then, along the semigroup, the following 
bounds hold 

J ht I Ck log ht P dfi < J ho log hodfi Vfc = 0, N, 

where ht = ft/ p and ft is the law of the process with generator L = —C. 
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B Ergodicity 

In this appendix we apply Markov chain techniques [^[^[5^ to prove ergodic- 
ity of the Markov process x(i) := {q{t), p(t), z{t)} given by^|). To simphfy the 
notation, we set all constants equal to 1. We will consider the ergodic properties 
of the SDE 

q = p 

p = -VqV{q)+r 
r = —p — r + W. 

We study either the case q or q T"^ and p, r S M'*; we will sometimes use 
the notation x = (q-,p , r) and L — —C is the generator of the process. 
Motivated by [2S], consider the following conditions: 

Lyapunov Condition : There exists a function G{x) : B?"^ — >■ [l,oo) such 
that G{x) — > cx) as ||a;|j — ?► oo and LG{x) < —aG{x) + d for some a, d > 0. 

Minorization condition : Let Pt{x, A) be the transition kernel of the Markov 
process x(t). There exist T > 0,7y > and a probability measure v, with 
i/(C"^) = and v{G) = 1 for some fixed compact set G in the phase space, 
such that 

Pt{x,A)> nv{A) \fA(^B{W^'^),x eG. 

We will use the following result, whose proof can be found in 

Theorem B.l. ; Minorization condition + Lyapunov condition ^ ergodicity. 

Assumption (★): Let Bs{y) £ R^'^ be the ball of radius s centered in y. For 
some fixed compact set G we have 

• Pt(x, A) has a density pt{x, y) which is continuous V(x, y) £ G x G, more 
precisely 

Ptix,A)^ f pt{x,y)dy VA e ^(R^'*) n 6(C), e C; 

J A 

• Vx £ G and V(5 > one can find a. t — t{6) such that 

Pt{x,Bs{x* j) > for some x* £ int{G),^x £ G. 

Lemma B.l. Assumption f+j Minorization Condition. 

Now let V{q) be any C^(T'*) potential, V{q) > —k for some positive constant 
k. Consider the function 

G{x) = C + f + + DV{q) + H{p, r), 

where B,G,D,H and C are positive constants to be chosen. We have that 

Gix) >c+%r+^\\rr-^\\pr- ^wrr Dk, m 
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so we need B > H , C > H and C > Dk. Moreover 

LGix) = D{VqV,p)-B{VgV,p)-H{VgV,r) + B{r,p)+H\ 
-C{p, t) - H\\pf - C||r||2 - H{p, r) + C 
H 
1 



< H\\rf + -||V,F||2 - H\\pf - C\\rf + C + H\\r\\\ 



where we have chosen B — D — C + H . On the other hand, since V(q) < K, 
G{x) > -^B\\pr ^CWrf aKB - a^Wrf - a^Wpf 

so imposing also 2H - C < -f (C + H), ~H < ~f (B + H) for some a > 0, 
the Lyapunov condition is satisfied. One possible choice is a = 1/4, B = 13/16, 
C = 5/8 and H=3/16. 

Notice that from what we have just proven it follows that VZ > 1 we have 

LG{xy <-aiG{xy +di, (49) 
for some suitable positive constants a; and di. In fact, 

dg^G{xy = lG{xy-^dq^G{x), 

dp^G{xy = lG{xy^'^dp^G{x), 

and 

d^Gix) = dr^ [lG{xy-^dr^G{x)] = l{l-l)G{xy-'^{dr^Gf +lG{xy-^dlfl{x). 
Furthermore, using p8)) . we obtain 

l{l - l)G{xy-''{dr,Gf < ciGixy-\ 

so that 

LG{xy < lG{xy-^LG{x) + ciG{xy-^. 

Hence, using what we have proven in the case ^ = 1, we get 

Consider now the case x(i) g M'' x M*^ x R'^. We introduce the Lyapunov 
function 

G{x) = C + ^\\qf + ^\\pr + ^\\rf + DV{q) 

+Eip, q) + F{q, r) + H{p, r) + M{VqV, p). 

Consequently, 

VqG = Aq + DVgV + Ep + Fr + MVV(g) • p, 

VpG = Bp + Eq + Hr + MV gV, 
VrG = Cr + Fq + Hp. 
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Thus, 



LG{x) = A{p,q) + D{W,V,p)+E\\p\\^ + F{p,r)-B{WgV,p) 

~E{V,V, q) - HiVqV, r) + B{p, r) + E{q, r) + H\\rf 
-C{p, r) - F{p, q) - H\\p\\^ + M{p, HessV{q) ■ p) 
-C\\rf - Fir, q) - H{p, r) + C - Af || V^^f + Af(r, VgV) 

From Assumption 12 . 1 f iii) it follows that there exist constants /3 and a such that 

ctIIqIP - /3||V,V^|p +00 as WqW^ +00. 

Hence, it follows that G satisfies the Lyapunov condition. Also in this case, one 
can prove that the Lyapunov condition holds for G(a;)', Z > 1, as well. 

As for Assumption (★), first of all let us notice that, since the operator 
dt + C is hypoelliptic, the transition probability has a density; because the SDE 
we are considering has time independent coefficients the density is C°° provided 
that V{q) € C°° |32| . Moreover, studying the control problem associated with 
dx — b{x)dt + adw, namely dX — h{X)dt + adU where U{t) is a smooth control, 
and using Strook-Varadhan support Theorem, we can prove that Pt{x,A) > 

OVx G R3^, t > and for any open A e R^'^. 

Consider now the set Qi = {g : M^d _^ measurable :| g{x) \< V{xy}. We 
have proven that the process x(t) has a unique invariant distribution p. Fur- 
thermore, there exist constants k — k{l) and A — A(Z), such that Vg E Qi 

\E^^g{x{t))~p{g)\<kG{xo)e-^\ t>0. (50) 
This completes the proof of the theorem. 
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